Simulation of Einstein-Podolsky-Rosen experiments in a local hidden variables model 

with limited efficiency and coherence. 



O 

o 

(N 



O 

m 

> 

o 
m 
o 

o 

Oh 

-I— > 

a 
cr 



W. A. Hofer 

Department of Physics and Astronomy, University College, Gower Street, London WC1E 6BT, UK 

We simulate correlation measurements of entangled photons numerically. The model employed is 
strictly local. The correlation is determined by its classical expression with two decisive difference: 
we sum up coincidences for each pair individually, and we include the effect of polarizer beam 
splitters. We analyze the effects of decoherence, detector efficiency and polarizer thresholds in 
detail. The Bell inequalities are violated in these simulations. The violation depends crucially on 
the threshold of the polarizer switches and can reach a value of 2.0 in the limiting case. Existing 
experiments can be fully accounted for by limited coherence and non-ideal detector switches. It 
seems thus safe to conclude that the Bell inequalities are no suitable criterium to decide on the 
nonlocality issue. 
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The Einstein-Podolsky-Rosen (EPR) problem has long 
occupied a central place in the understanding of quantum 
mechanics 0], Bell's inequalities in conjunction with cor- 
relation measurements seemed to prove that reality in mi- 
crophysics is manifestly nonlocal Furthermore, the 

experimental evidence seems to contradict even the no- 
tion of an independent reality M . Both of these features, 
if true, are highly problematic. The former, because no 
field propagating with a velocity exceeding c has ever 
been observed. The latter, because without an indepen- 
dent reality there is no guarantee that theoretical models 
can at all be contradicted. And without the possibility 
of contradiction scientific progress follows no clear rules. 

For these reasons the EPR problem is far more impor- 
tant than the experiments alone indicate. Consequently, 
a large amount of work has been devoted to this field. 
Two years ago, the standard reference on EPR - the book 
by Afriat and Selleri [|| - more or less highlighted the 
dilemma. But in the same year Deutsch and Hayden 
could show, by an analysis of the information flow 
in such an experiment, that there is in fact no nonlocal 
connection between the two measuring devices. All in- 
formation about the two angles of polarization, <f>i and 
4>2 , is stored locally. Even though this information can- 
not locally be accessed. It is probably due to this new 
field of research, quantum information theory, that the 
problem is even more important today than it was ten 
years ago. Consequently, a number of papers in the last 
two years have analyzed the paradox from different an- 
gles, and the analysis brought two features into focus: (i) 
The validity and significance of Bell's inequalities 
and (ii) the relevance of a photon's phase for the corre- 
lations ]lC|-|T^] . From the viewpoint of information flow 
a violation of Bell's inequalities is no proof of nonlocal- 
ity [0. From a formal point of view it seems that the 
standard inequalities cannot be derived without violat- 
ing established notions about the measurement process. 
For a detailed discussion see Sica or Adenier |||)]]. The 
notion of a phase seemed initially problematic because 



the phase e.g. of a wavefunction cannot be related to 
physical properties of a photon. But as shown later, the 
phase indicates the phase of a photon's electromagnetic 
field jll]] . And it could be established that the existence 
of a phase connection between the two points of mea- 
surement, a connection which arises at the process of 
emission from a common source, is sufficient to explain 
correlations between two measurements in space-like sep- 
aration. It was also emphasized that measurements can- 
not in general be factorized without loosing the linearity 
of the fields between the two polarizers. 

In this Letter we pursue a different strategy. We per- 
form numerical simulations of actual experiments. We 
sought to include the features of the experimental sit- 
uation as far as possible. For this reason we shall give 
not only the results of ideal measurements, but also mea- 
surements with limited efficiency, limited coherence, and 
under the condition of dead angles of our polarizers. It 
will be seen that all these effects have a bearing on the 
actual data. 

Setup. - The setup of the experiment is shown in Fig. 
|l| (a). A source between two polarizers emits a pair of 
photons along the z axis. Two polarizers, at the posi- 
tions L\ and L2, respectively, measure the angle of po- 
larization. The angle of polarizer one is varied by a half 
cycle, 7T, during the experiment. At every position of 
the polarizer a set of 1000 photon pairs is emitted and 
measured. The switch of the polarizer is shown in Fig. 
[I] (b). If cos 2 (0i — a) is larger than 0.5 + As, then the 
event is recorded as a transmission (+). If it is less than 
0.5 — As, it is recorded as an adsorption (-). No ac- 
tion is taken for values between these two boundaries. 
The same switch is applied to both measurement devices. 
The threshold As formalizes three separate features of 
the experiments: (i) The transmission characteristics of 
the polarizer beam splitters; (ii) the threshold of the pho- 
todetectors; and (hi) the electronic evaluation of events, 
since double counts (+ and —) are excluded. 

A single simulation run starts with the initialization of 
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the random number generator |14| . The generator is ini- 
tialized only once, at the beginning of a full simulation 
cycle. The random number is mapped onto the initial 
phase from to 2tt of the photon pair. Simulations are 
generally made with a phase difference of tt/2 between 
the angles of polarization of photon one and photon two. 
After covering the distance to polarizer one and two the 
photons are measured. We assumed, without lack of gen- 
erality, that both distances are integer multiples of the 
wavelength. After a single pair has been measured, we 
record the coincidences (++,+-,-+,-). The procedure is 
repeated for all 1000 pairs, then the polarization angle 
of device one is changed by 7r/100. A run ends, when all 
1000 pairs at the final position of polarizer one have been 
measured (it) . In all figures we only plot the coincidences 

We accounted for limited efficiency and decoherence in 
the following way. Limited efficiency means that not all 
pairs emitted are actually measured. In this case we sim- 
ply did not evaluate all pairs, depending on the efficiency 
of the setup. 50 % efficiency, for example, means that 
only every second pair is actually recorded. To simulate 
decoherence we created an independent random input for 
a certain fraction of a half cycle of tt. 100 % decoherence 
here means that half a wavelength of the photon's opti- 
cal path is random. This translates into a polarization 
angle random in the interval [0,7r]. Both effects reduce 
the maximum of the output measured, but it will be seen 
presently that they have very different effects. 

Ideal measurements. - Initially we simulated an ideal 
measurement. The efficiency in this simulation is 100%, 
the fields of the two photons are fully coherent through- 
out the distance between the two polarizers, and the ex- 
perimental devices are supposed to have ideal character- 
istics. The result of this simulation is shown in Fig. |^. 
We did two separate simulations, one with a polarization 
difference between the two photons of (full squares), 
the other with tt/2 (full circles). It can be seen that 
neither of the simulations comes close to the theoretical 
prediction of a sin 2 ((3 — a). Instead, the curves repre- 
senting ideal measurements would be of angular shape. 
However, the maximum of the correlation (No/2, where 
Nq is the number of pairs) and the minimum (zero) are 
exactly obtained in the extreme cases. 

It should be noted that the results given in these plots 
reflect the "classical" formulation of a coincidence, given 
by the equation EH]: 

P(a,/3) = J dXcos 2 (X- a)cos 2 (X- (3) (1) 

with two decisive differences: first, the summation is 
performed over single pairs, as in the actual experiments, 
rather than over the two polarizers separately. The latter 
procedure, formalized in the given integral, includes not 
only photons of one pair, but also sums up contributions 
of different pairs. And second, we accounted for the digi- 
tal output of the polarization devices. The integral is no 



suitable representation of the digital results in current 
experiments. If, for example, X — (pi and <p2 = <pi + tt/2, 
then for a = (3 — the integral yields No/8. But the 
actual count, under the condition that cos 2 (pi > 0.5 and 
cos 2 (</>i + tt/2) > 0.5 is zero. The digitalization, neces- 
sary to obtain formal agreement with spin measurements 
in quantum mechanics, is usually achieved by means of 
a polarizer beam splitter [|J. 

Dead angles. - In our simulations we find that the 
curves obtained are not very sensitive to the threshold of 
the polarizer switches. We have performed simulations 
where As was varied from 0.00 to 0.20. Apart from a re- 
duction of the absolute yield the increase of the threshold 
only affects the width of the minimum at the ultimate an- 
gles. This effect is equal to a retardation of the onset of 
the correlation function at its minimum position. The 
threshold therefore does not change the functional form 
of the correlations. 

Efficiency. - The detection of photons is one of the 
problems experimenters are still confronted with. The 
efficiency is in fact so low (less than 10% 0), that the 
correlations found in Aspect's measurements p] were dis- 
puted on the grounds of an "efficiency loophole". In 
our simulations such a conclusion would not be justi- 
fied. Even though the shape of the curve changes some- 
what and the statistical spread is dramatically increased 
in the low efficiency range, the maxima and minima are 
still clearly distinguishable. The minimum, moreover, re- 
mains zero. We included a dead angle of detectors by a 
threshold of 0.05. The same threshold will be used in 
subsequent simulations. The coincidence rates due to 
detector efficiency varying from 50% to 10% are shown 
in Fig. |3|. From this result we conclude that efficiency 
is not the decisive issue to estimate the relevance of an 
experiment. 

Decoherence. - The polarizers in current measurements 
are more than 400 m apart Furthermore, there is 
no vibration damping or cooling to very low tempera- 
tures involved in such a measurement. This feature of 
the measurements is bound to cause random motion of 
system components. From surface science the range of 
motion without damping can be estimated, it should be 
for an isolated surface no less than a few nanometer or 
more than one percent of the photon's wavelength. Con- 
sidering that we deal with three coupled components and 
optical paths in between it seems safe to increase this es- 
timate by one order of magnitude. In this case we have 
to include random motion of our system in the range of 
about 5-10% of the wavelength. This translates, in our 
simulation, into a rate of decoherence of 10-20% (100% 
means that half a wavelength of the photon's optical path 
is random). 

Simulation with a decoherence rate of 10%, 50% and 
100% are shown in Fig. [|. The interesting feature of 
decoherence is that it renders the resulting distribution 
more sinusoidal than the correlations of an ideal mea- 
surement. The fully decoherent simulation proves that 
correlations are independent of the setting as required, 
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but it also shows the noise due to a random distribution 
of the initial phase of the coupled system. In practice 
all effects analyzed will to a greater or lesser extent be 
present in any single measurement. 

Bell violations. - Finally we demonstrate the influ- 
ence of the polarizer threshold on the violation of Bell's 
inequalities. To this end we simulate the counts at four 
selected angular positions of the polarizers a and j3 (0°, 
45°, 22.5°, 67.5°). These positions yield the maximum 
violation of Bell's inequalities in the standard framework. 
We performed the simulations for varying threshold val- 
ues from 0.0 (no threshold) to 0.2 (nearly half the photons 
remain undetected). For every setting we performed 10 
separate runs, each with 10000 pairs of photons, the ef- 
ficiency of the detectors was assumed to be 100%. Fig. 
^ gives the result of our simulation. The violation (com- 
puted according to the version of Clauser et al. (CHSH) 
p6|) increases with increasing threshold. Furthermore, 
the limit of violation is close to 2.0 (CHSH value of 3.90) 
in the final setting. The actual threshold in the experi- 
ments can be estimated from the visibility of the corre- 
lation function. For decoherence rates of 10 - 20 %, as 
in the experiments, we can only obtain agreement with 
reported values (97% visibility Q), if we increase the 
threshold As to more than 0.10. Our simulations indi- 
cate that the violation depends crucially on the thresh- 
old. The maximum violation can reach a value of as much 
as 2.0 (the limiting case). It seems thus safe to conclude 
that the Bell inequalities are no suitable criterium to de- 
cide on the nonlocality issue. 

Summary - In summary we have presented a numerical 
simulation of EPR experiments under the assumption of 
strict locality and analyzed the effects of polarizer thresh- 
olds, limited efficiency, and decoherence in detail. Wc 
could show that Bell's inequalities are violated in these 
simulation, and that the violation depends crucially on 
the threshold of the polarizer switches. The limit of vio- 
lation in this model is about 2.0 (CHSH value of 4.0). We 
found that existing measurements can be fully accounted 
for by limited coherence and non-ideal polarizer switches. 
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FIG. 1. One dimensional model of EPR type experiments, 
(a) The measuring devices (1) and (2) are in opposite direc- 
tions from the photon source S. The polarizers are set to the 
angles a and /3, respectively, (b) The switches at both sta- 
tions measure the polarization and, depending on the angle, 
either transmit 1 or to the computer. Note that the dead 
angle of the polarizer is simulated by a threshold As of the 
switch. 
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FIG. 2. Ideal measurement of correlations (N ++ ). Two 
simulations were performed with a difference of (full 
squares) and 7r/2 (full circles) of the polarization angle at 
the origin. Neither of these curves is equal to the theoret- 
ical prediction, instead we obtain an angular shape for the 
correlations. 
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FIG. 3. Dependence of coincidence rates on the efficiency 
of detection. The shape of the distribution is similar to the 
ideal distribution, but the statistical spread is considerably 
larger. We include the sin 2 (a) function for reasons of compar- 
ison. This function has not been actually fitted to simulated 
data. 
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FIG. 5. Violation of Bell's inequalities depending on the 
threshold of the polarizer switches, (a) The inequality is vio- 
lated in all cases where the threshold is not zero (full circles). 
We obtain a maximum CHSH value of 3.9 (threshold 0.2). 
The experiments of Weihs et al. indicate a threshold of 0.1 
(empty circle), (b) Simulation of EPR experiments with a de- 
coherence of 10% and a threshold of 0.1. It can be seen that 
the distribution differs only insignificantly from the square of 



FIG. 4. Dependence of coincidence rates on the decoher- 
ence of photon beams. Due to decoherence the distribution 
becomes more sinusoidal. In the limit of full decoherence we 
obtain uncorrelated measurements. 
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